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1 . 0  Introduction 

The  United  States  Army  has  recognized,  especially  as  a  result  of  tactical 
communications  experience  in  Vietnam,  the  need  for  a  realistic  model  for 
wideband  UHF  radiowave,  propagation  within  different  types  of  forest  [1], 
Although  the  early  dielectric  slab  models  [2,3,4]  correctly  explain  the 
up -over -and -down  lateral  waves  propagating  at  the  air- forest  interface  over 
long  VHF  paths,  they  characterize  the  forest  only  in  terms  of  an  effective 
dielectric  constant  which  is  postulated  rather  than  derived  from  the 
biophysical  characteristics  of  the  vegetation.  Further,  the  dielectric  slab 
models  are  not  applicable  at  UHF  and  above  where  incoherent  scatter  dominates 
and  wideband  signals  suffer  appreciable  delay  spread. 

More  recently,  CvberCom.  with  the  support  of  US  Army  CECOM  and  the  co¬ 
operation  of  SRI  International,  has  undertaken  the  development  of  a  realistic 
stochastic  propagation  model  for  characterizing  the  wideband  forest  channel. 
This  program,  both  theoretical  and  experimental,  addresses  several 
interrelated  study  areas:  the  biophysical  characterization  of  the  forest; 
electromagnetic  interaction  of  the  propagating  radiowave  with  the  scat.terers 
using  theories  of  discrete  multiple  scattering  and/or  transport  theory;  and 
communications  channel  characterization  in  terms  of  transmission  loss  and 
delay-Doppler  spread.  The  early  CvberC'om  studies  were  restricted  to 
trunk- dominant  coniferous  forests  [5,6,71;  more  recently,  canopy- dominant 
deciduous  forests  have  been  considered  [8], 

The  acquisition  of  experimental  data  supporting  the  development  of  a 
propagation  model  for  canopy- dominant  deciduous  forests  began  in  August,  1987 
within  an  experimental,  red  maple  forest  located  near  Coventry,  Connecticut. 
Anticipating  that  model  development  and  validation  might  require  a  detailed 
biophysical  description  of  the  forest,  this  site  was  selected  because  it.  had 
been  the  subject  of  a  recent,  detailed,  quantitative,  biophysical  study 
conducted  by  its  owner  and  manager,  the  University  of  Connecticut. 

The  experiments  at  Coventry  were  first  conducted  during  the  summer's  end 
of  1  987  (24  August  to  4  September),  and  then  repeated  in  November  after  the 
autumnal  fall  of  the  leaves  in  order  to  assess  the  seasonal  effect  of  the 
foliage.  Although  data  from  the  first  set  of  measurements  had  been  processed 
and  analyzed  [9],  data  from  the  second  set:  had  not  .  A  selected  sub-set  of 
these  data  is  analyzed  in  this  report  to  assess  the  effects  .(  foliage,  rain, 
and  snow  on  delay-spread  in  canopy -dominated  deciduous  forests . 
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2 . 0  Background 
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From  the  communications  viewpoint,  radiowave  transmission  through,  forests 
may  be  modeled  by  a  randomly  time -variant  linear  channel  whose  input  and 
output  signals  are  related  by  the  integral  transform  [10] 

y(t)  -  Jx(t-Og(t,  ?)d£  (2-1) 

where  x(t)  is  the  transmitted  (input)  signal,  y(t)  is  the  received  (output) 
signal,  and  g(t,|)  is  the  response  cf  the  channel  at  time  t  to  a  unit  impulse 
applied  at  time  t-£.  Thus,  the  stochastic  behavior  of  the.  channel  is 
manifested  through  g(t,£)  in  two  dimensions  via  the  propagation  delay  variable 
£  which  characterizes  the  spatial  configuration  of  the  forest  scatterers  and 
via  the  time  variable  t  which  characterizes  their  (possibly)  time-variant 
configuration  Since  x(c)  and  y(t)  are  scalars,  the  vector  coupling  of  the 
forest- scattered  electromagnetic  fields  to  the  transmitting  and  receiving 
antennas  is  assumed  to  be  implicitly  embedded  within  g(t,£). 

In  order  to  characterize  the  stochastic  behavior  of  the  forest  channel 
completely,  multi -dimensional  probability  density  functions  of  g(t,£)  are 
required.  However,  because  of  multiple  scattering  within  the  forest,  the 
received  signal  may  be  considered  to  arise  from  the  superposition  of  a  large 
number  of  independent  contributions  having  the  same  relative  delay  As  a 
consequence  of  the  central  limit:  theorem  it  may  be  concluded  that  y(t.)  is  a 
gaussian  random  process  and,  in  view  of  Equation  (2-1),  so  too  is  g(t,£).  The 
principal  utility  of  this  hypothesis  lies  in  the  consequence  that  all  orders 
of  the  multi -dimensional  probability  density  function  describing  g(t , £)  can  be 
expressed  in  terms  of  the  correlation  function 


If  the  channel 
s  t  a  t  i  o  n  a  r  y  I  n  the 
d i  f  fere 1 i 1:  de  1  ays 


R  (t  t  •£  ,£  )  -  <g(t  J  )g*(t 

K  1  Z  1  Z  1  2 


f  _  )  > 


transfer  function  g(t,F)  can  be  considered 
time  variable  t  and  the  scattered  signals  corres 
f  uncorrelated,  then  the  channel  correlation 


(2-2) 

wi de -  sense 
ponding  to 
f uno t i on 


simplifies  to 


,  ( t ,  f  )  g  ( t:  t  t  ,  t) )  >  -  Q  (  r  ,  £  )  6  (  f  ■  t] ) 


( 2  -  3 ) 


where  6  ( • )  i.s  the  Dirac  delta  function  and  Q(  r,()  is  t  tie  co- variance  f  unc  t  ion 
of  those  channel  fluctuations  having  propagation  delays  in  the  interval 
( t  ,  f'+d.f  .)  .  The  function  Q(0,f)  Is  variously  cal  led  delay  power  spent  turn  or  the 
del  a  y  -  s  p  r  e  a  <  i.  t  u  n  c  t  ion. 
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In  the  absence  of  relative  forest  motion  due  to  wind  and/or  terminal 
movement,  the  forest  channel  [and  g(t,£)j  may  be  considered  time  -  invariant . 
The  ensemble -averaged  channel  correlation  function  then  simplifies  to 

<g(t,Og  (t.*7)>  =  Q (£)S(£-n)  (2-4) 

whe  re 

Q(f)  -  Q(0,€)  (2-d) 

is  the  delay-spread  function. 


A  physically  appealing  interpretation  of  the  delay-spread  function  can  be 
derived  by  considering  the  received  power  when  an  unmodulated  sinusoid  is 
transmitted.  In  this  case,  the  transmitted  signal  is 


x(t)  =  exp{j2?rf  t) 

O 

so  that,  according  to  Equation  (2-1),  the  received  signal  is 

y(t)  -  I"  expt  j2mfQ(t-£)  )g(t,£)d£ 
and  the  received  power  is 

<P  >  =  <yZ(t)>  -  ff  exp  { j  2  7r  f  ( t:  -  ,f  )  )  exp  (  -  j  2.n  f  ( t.  -  £  )) 
r  J  J  o  1  o  2 

<g(t:-?,  )g*(t.?2)>dCx<if2 


f  exn{-i2flf  (f  -f  )  )Q(f  )S(£  -£  )df  df 
J  J  ‘  ■'  o  i  ‘  2  1  i  z  ‘  1  ‘ 


(2-6) 


(2-7) 


(2-8) 

(?.-<>) 
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-  f  Q(£)d£  (2-10) 

j 

Since  £  denotes  the  propagation  delay,  it  is  apparent  from  Equation  (2-10) 
that  Q(£)  describes  the  ensemble  -  aver  a  ged  oist  ri but  ion  of  the  Increment  a l 
power  received  within  n  delay  interval  d£ .  According  to  the  forest  propa¬ 
gation  model  develop-vl  in  reference  [  /  J  ,  the  theoretical,  en.semb  -averaged 
de lay  -  spread  of  a  homogeneous,  t rank - domi nant  forest  is  of  the  loim 
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signal  is  impulsive,  i.e. 


R  (r)  =  <x(t)x(t-r)>  3  <S(r)  (2-12) 

X. 

Cross-correlation  of  the  received  signal  with  x(t-r)  yields 

<X  (t-r  )y(t)>  -  J  <x(t-r)x(t-^)>g(t,Od^  -  J  R^r -f  )  g(  t ,  £  )  d£  -  g(t  ,rj  (2-13) 
so  that  to  the  extent  that  R  (r)  is  approximated  by  a  Dirac  delta  function, 

X 

the  cross  -  correlated  output  g(t,r)  approximates  the  input  delay -spread 
function.  The  delay- spread  function  Q(»)  is  approximated  by 

P(r)  =  I g( t , r )  |  2  (2-14) 

This  measured  counterpart,  distinguished  here  by  the  name  Power  Time -Variant 
Impulse  Response  [PTVIR],  is  shown  in  Figure  2-2.  Its  ragged,  stochastic 
structure  can  be  attributed  to  phase  interference  between  forts t- scattered 
radiowaves  having  the  same  propagav ion  delay  r  and  so  suggests  that  the  power 
P(r)  is  an  exponentially-distributed  random  variable. 

For  time -variant  ergodic  media,  t  iine  -  averaging  can  be  used  to  smooth 
successive  PTVIR  measurements  for  comparison  with  their  theoretical 
ensemble- averaged  counterparts.  However,  except  lot  terminal  and/or  wind  • 
induced  canopy  motion,  the  forest  channel  is  essentially  t  i  me  i  nva  r  i  a  nt  .so 
that  t  ime  -  averaging  plays  no  role.  Although  the  PTVIR  might  st  ill  lie  .smoothed 
by  curve  -  fitting,  t  he  performance  of  digital  radio  communion  t.  ion  systems 
performance  is  not  particularly  sensitive  to  the  no  of  the  [’TV  I R  [l.1!]  but 

rather  only  to  its  width  or  "delay  spread".  Theiefore,  a  sut  iirieut  (and 
certainly  simpler)  comparison  between  t  iieorv  and  experiment  can  he  eifeetod  in 
terms  of  the  delay  spread. 
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3 . 0  Experiment. 


m 


3 . 1  Forest  Site  Description 

The  Coventry  site,  a  flat  floodplain  with  a  30-year  old  stand  of  red 
maple  [Acer  rubrum]  with  occasional  white  pine  [Pinus  strobus]  and  trembling 
aspen  [Populus  tremul.oides  ]  at  the  edges,  is  part  of  the  University  of 
Connecticut's  forest  research  station.  A  photograph  and.  plan  of  the  Coventry 
site  can  be  found  in  Figures  3-1  and  3-2,  The  measured  propagation  paths  are 
identified  on  the  site-  plan. 

A  regular  grid  of  20  circular  plots  of  20-fcot  radius  was  laid  out  with 
100- foot  separations.  The  number  of  trees  and  their  diameters  within  each 
plot  were  recorded  and  used  to  calculate  the  tree- trunk  number  density  and 
mean  tree -trunk  diameter  which  were  then  used  r.c  construct:  the  contour  plots 
shown  in  Figures  3-3  and  3-4  and  the  trunk-diameter  histogram  of  Figure  3-5. 
The  stem  density  path  profiles  shown  in  Figure  3-6  were  derived  from  the 
tree- trunk  number  density  contour  plots.  Using  a  platform  elevator  provided 
by  the  University  of  Connecticut,  the  mean  tree  height  was  estimated  to  be  4? 
feet  (14.3  m)  with  no  more  than  about  5  percent  of  the  trees  exceeding  61  feet 
(18.6  m) ;  the  live  canopy  began  at  a  mean  height  above  ground  of  about  35  feet 
(10.7  m) . 


The  deciduous  canopy  characteristics,  recently  measured  by  Professor 
David  Miller  and  his  students  [13,14,15],  are  summarized  in  in  Figures  3-7, 
3-8  and  3-9  which  show,  respectively,  leaf -area- index  (LAI)  and  histograms  of 
leaf  inclination  angle  and  leaf  azimuth  angle.  The  leaf  number  density  p  and 
the  fractional  volume  occupied  by  the  foliage  f  can  be  estimated  from  the 
re lat ions 


LAI 

p  \  A*h 


(3-1) 


f  -  p  *  V 

V  1  1 


LAI  •  v 
h 


(3-7) 


where  LAI  denotes  the  leaf  -  area  -  index ,  A  is  the  average  area  per  leaf,  t  is 
the  average  leaf  thickness , and  h  is  the  canopy  thickness.  Typically,  for  red 
maple,  A  is  25  cm  and  t.  is  0  2  mm.  Thus,  for  ■  measured  LAI  of:  4.5  (refer  to 
Figure  3-7)  and  canopy  thickness  of  3  6  in  (see  above.)  ,  the  corresponding  leaf 
number  density  is  500  leaves  per  cubic  meter  and  :  .e  leaf  fractional  volume  is 
0.023  percent. 
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igure  3-4:  Coventry  Average  Trunk  Diameter  Contours  (1SB7) 
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igure  3-6:  Stem  Density  Path  Profiles  (Coventry,  1987} 
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Figure  3-8: 


Coventry  Leaf  Inclination  Angle  Histogram 
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3 . 2  Channe 1  Probe 

The  experimental  radiowave  propagation  data  were  acquired  using  the  US 
Army's  Wideband  Propagation  Measurement  System  (WPMS) .  The  WPMS ,  designed  and 
built  for  the  US  army  by  SRI  International  [16,17],  was  developed  especially 
for  measuring  ground-to-ground  communication  channel  characteristics  (trans¬ 
mission  loss  and  delay-Doppler  spread)  in  the  200-2000  MHz  band.  This 
computer-controlled,  mobile,  radio- channel  probe  uses  a  direct  -  sequence , 
pseudo  -  randomly  modulated,  wideband  waveform  and  sliding  correlator 
architecture.  Its  capabilities  are  summarized  in  Table  3-1. 

The  essential  elements  of  WPMS  operation  can  be  described  with  the  aid  of 
the  diagram  shown  in  Figure  3-10.  The  transmitter  sends  out  a  repetitive 
pseudo  -  random  waveform,  X^,(t),  with  a  "chip”  interval,  At,  and  a  code  period, 
T.  The  in-phase  (I)  and  quadrature -phase  (Q)  components  of  the  received 
signal  are  then  cross -correlated  against  a  replica  of  the  transmitted 
pseudo- random  waveform  in  each  of  four  parallel  sub -channels  successively 
offset  from  each  other  by  one-quarter  chip  (to  improve  TVIR  resolution).  The 
output  of  each  correlator  is  integrated  for  T  seconds  and  then  sampled.  The 
receiver  pseudo  -  random  code  generator  is  Chen  delayed  (slipped)  one  chip  and 
the  integrate-and-sample  operation  repeated.  The  complete  cycle  is  repeated  N 
times,  where  N  is  the  number  of  chips  in  the  pseudo  -  random  code  period.  Thus, 
each  complex  receiver  sub-channel  produces  one  output  every  T  •+•  At  seconds. 
To  obtain  a  /alid  measurement  of  the  entire  delay-spread  function,  the  channel 
must,  remain  essentially  time  -  invariant  for  N(T-t-At)  seconds.  The  resulting 
output:  is  stored  and  displayed  as  a  time-varying  impulse  response  (TVIR)  of 
length  MAC  and  resolution  At. 
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Table  3-1:  WPMS  Capabilities 


General 

Carrier  frequency  range 

Delay- spread  range 

Delay-spread  resolution 

Doppler- spread  range 

Doppler- spread  resolution 

TVIR  amplitude  resolution 

TVIR  multipath  amplitude  resolution 

Measurable  path  loss 


Antenna  3 

Omnidirectional 

Directional 

Transmitter  polarization 
Receiver  polarization 
Azimuthal  range 
Height  range 


Transmitters 

Number 

frequency  range 
XMTR  no.  1 
XMTR  no.  2 
Power  output 
Modulation 

DS  code  length  (chips) 

DS  code  race 

Bandwidth  (null- to -null) 

Cont  rol 


Receivers 

Number 

Frequency  range 
Noise  figure 
Demodulation 


DS  code  length  (chips) 

DS  code  rate 

Bandwidth  (null-to-nui 1 ) 
Control 

Input  signal  range 
Instantaneous  dynamic  range 
AGO  (computer  controlled) 


200-2000  MHz 
1-20  ^-sec 

1  nanos  ic 
15-240  Hz 

2  Hz  (max) 
0.1  dB 

-20  dB 

155  dB  (max) 


Biconical 

Log-Periodic  Dipole 
V,  H,  RCP 
V,  H.  RCP,  LCP 
360  degrees 
5-65  feet 


2 

200-1050  MHz 
700-2000  MHz 
100  W  (max) 

Direct- Sequence  BPS: 
Continuous  Wave  (CW) 
255,  511,  1023  or  2047 
50,  125  or  250  MHz 
100,  250  or  500  MHz 
Local  or  Remote  (RCVR) 


2 

200-2000  MHz 
9  dB 

Co he rent 

Direct-Sequence  BPSK 
Continuous  Wave  (CW) 
255,  511,  1023  or  2047 
50,  125  or  2.50  MHz 
100,  2 50  or  500  MHz. 

Local  compute i 
-95  -  0  dBm 
50  dB 

L  dB  increment:.*.: 


T/4 


8  192 
I 


Rico  rcis 


3 . 3  Data  Transfer 

Unfortunately,  the  HP- 1000  computer  is  now  obsolete  and  the  computer 
facilities  available  to  CvberCom  [IBM  370  and  HP-9000]  were  unable  to  read  the 
binary  data  tapes.  Although  a  number-builder  might  have  been  constructed, 
CvberCom  elected  to  use  an  HP- 1000  still  available  at  CECOM  in  Fort  Monmouth, 
New  Jersey  to  read  the  tapes  and  to  translate  the  data  into  ASCII  format  for 
compatibility  with  the  IBM  370  [and  most  other  computers].  Although  the 
HF-1000  software  used  by  SRI  International  for  reading  the  binary  tapes  was 
made  available  to  CvberCom.  these  programs  proved  unwieldy  to  use  and  provided 
superfluous  data  detail.  A  simpler  program,  developed  by  the  National  Bureau 
of  Standards  and  acquired  from  them,  was  used  instead  [see  Appendix  A], 

Earlier,  CvberCom  had  developed  some  of  the  software  that  would  be  used 
for  propagation  data  analysis.  These  programs,  however,  ran  only  under 
HP-BASIC  on  CvberCom' s  HP-9816  computer.  Unfortunately,  this  computer  does 
not  interface  with  the  nine-track  tape  drives  required  to  read  the  HP- 1000 
ASCII  data  tapes.  However  these  tape  drives  were  available  locally  at.  the 
George  Washington  University  in  Washington,  D.C.  where  CvberCom  maintains  a 
computer  account.  It  was  decided  then  to  transfer  the  ASCII,  data  from  tape  to 
disk  using  GWU's  IBM  370  computer  and  then  to  download  the  files  to  CvberCom' s 
HP-9816  using  telephone  lines. 

Using  IBM's  Job  Control  Language  (JCL)  the  automated  utility  program 
TRNSFR  [see  Appendix  B]  was  developed  for  GWU's  IBM  370  to  read  the  HP- 1000' s 
ASCII  data  tapes  and  re-write  a  reduced  data  set  [16  PTVIRs  from  the  first 
sub-channel  for  each  of  two  receivers]  on  disk  according  to  the  format  shown 
in  Figure  3-12.  The  reduced  data  sets  are  then  down- loaded  to  CvberCom  using 
cur  HP-9816  computer  and  HP-supplied  ASCII  Terminal  Emulator. 

At  Cyber Com  the  down- loaded  ASCII  data  are  temporarily  stored  on  3 . 5- inch 
floppy  disks.  The  limited  capacity  of  these  media  <,2^6  kilobytes)  Led 
CyberCom  to  develop  'he  software  utility  BCKUP  [see  Appendix  C]  for  re -storing 
the  ASCII  data  on  150-foot,  1/4- inch  wide,  16- track  magnetic  tape,  HP -88140SC 
data  cartridges  (16.7  megabytes/cartridge) .  Headers  identifying  salient  ex¬ 
perimental  parameters  are  prefixed  to  each  file  and  summarized  for  convenient 
access  in  a  Measurement  File  Catalog.  Another  program,  KELVIN  [.see  Appendix 
D),  was  developed  to  read  the  ASCII  data  from  the  tape  cartridges,  re- format 
it  for  consistency  with  prtr.  .ou;  ly  developed  analysis  so  I L  ware ,  and  re -store 
it  in  binary  (BOAT)  format  .  The  temporary  IibAT  fi  ie.s  are  used  in  data 
analysis  because  they  have  much  higher  transfer  rates  than  the  ASCII  files. 
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GWU  Reduced  Data  Set  File  Format 


© 


© 


© 


w 


Q 


© 


» 


The  time  to  transfer  the  data  of  a  single  episode  from  an  archived  file 
(HP-1000  BOAT  on  9 -track  tape)  to  a  reduced  data  set  (ASCII  data  on  1/4- inch 
tape  cartridge)  is  not  negligible  -  approximately  75  minutes  per  episode:  40 
minutes  are  required  by  the  HP-1000  to  convert  the  data  from  BOAT  to  ASCII,  5 
minutes  are  required  by  the  IBM-370  to  read  the  ASCII  data  and  re-write  the 
reduced  data  set  on  disk,  25  minutes  are  required  to  transfer  the  reduced  data 
set  to  CyberCom ,  and  5  minutes  are  required  to  re -store  the  data  on  tape 
cartridge . 

3 . 4  Data  Selection 

A  summary  overview  of  the  forest  experiments  can  be  found  in  SRI's  Final 
Report  [9],  CyberCom  reviewed  these  data  and  selected  a  subset  of  ten  (10) 
computer  tapes  representing  nearly  one -hundred  experiments  conducted  over 
different  distances  with  different  frequencies,  polarizations  and  antenna 
heights . 

3  .  5  Data  Aralys is 

The  quantitative  characterization  of  the  measured  PTVIR  P(r  ) 

complicated  by  the  stochastic  character  of  the  multipath  and  by  the 
contaminating  noise  (both  thermal  and  code)  contributed  by  the  radio 
equipment.  The  problem  is  illustrated  in  the  profile  of  Figure  3-13:  region 
A  is  the  noise -dominated  ''precursor"  region  attributable  primarily  to  thermal 
noise  in  the  receiver  front-end,  and  code  noise  due  to  imbalance  in  the 

correlators'  region  B  is  the  multipath-dominated  signal  region  where  noise  is 
negligible;  and  region  C  s  the  noise -dominated  "tail"  region  where  the 
multipath  signal  disa>  tears  into  the  noise. 

Because  the  received  signal  power  associated  with  any  particular 
multipath  delay  is  contributed  by  scattering  from  a  large  number  of  trees, 
P(r  )  is  an  exponent i ally- distributed  random  variable  (the  received  voltage 
then  being  Rayleigh- distributed) .  Fluctuations  in  F(r  )  can,  for  wide -sense 

k ' 

stationary,  time -variant  media,  be.  reduced  by  profile  averaging;  for  the 

time  -  invariant:  forest  channel  measured  on  a  single  pat  i  ,  profile  averaging 

serves  only  to  smooth  the  noise-dominated  regions  A  and  C  which  are  of  little 
intrinsic  interest,  «nd  does  not  provide  an  ensemble  average. 

Al  though  time- invariance  simplifies  data  acquisition  it  eompl  i  cat«  a  data 
ana  lysis  because  the  errodic  hypothesis  can  not.  be  invoice  tc  mouth 
measurements  by  t ime -averaging.  Fortunately ,  however,  radio  common l ca t  ion 
system  performance  is  not.  parti cu  arly  sen. si!  i ve  to  the  shape  >■(.  the  PTVIR 
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[12]  but  rather  only  to  Its  width  or  "spread."  Perhaps  the  most  universally 
accepted  and  tractable  measure  of  delay  spread  is  the  radius -of  gyration 
(square-root  of  the  sec.ond-centrai-moment)  defined  by 


T 

rog 

v  rn  m  -mm 

2  0  11 

/  ra 

0 

(3-3) 

where  the  moments  are  defined  by 

rT 

K 

m  = 
o 

Q(£)d£  - 

I  P(rk) 

(3-4) 

J  0 

0 

rT 

K 

m  ~ 

l 

£Q(Od£  - 

Y  t  P(r  ) 

L‘  k  ^  k7 

(3-5) 

’  0 

0 

.T 

K 

ttl2  "  j 

e2q - 

1 

(3-6) 

0 

0 

with  the  integrals  being  used  with  continuous  analytic  expressions  for  the 
delay-spread  function  Q(£)  and  the  summations  with  sampled-data  representa¬ 
tions  for  the  PTVIR  PCr^).  The  mean  delay  (center-cf -gravity)  of  the  PTVIR  is 
given  by 


avg 


m  /m 
1  n 


(3-7) 


The  "leading  edge"  of  the  PTVIR  where  the  received  power  first  becomes 
distinguishable  from  the  background  noise  of  the  receiver  will  be  denoted  by 


T  Excess  delay,  the  difference  between  T  and 

m  i.  a  av  g 


nun  ' 


1  S 


T  -  T  -  T 

exc  a  v  g  min 


(3-8) 


Still  another  measure  of  delay  spread  and  better  suited  to  the  non- 
symmetric  delay -spread  functions  measured  in  forests  is 


T  -  T  )  T 

spr  exc  r  o  g 


(3-9) 


t.omputer  program  MELVIN,  developed  to  calculate  t  he  moniciii  s  (in  }  is 

k 

listed  in  Appendix  K.  The  algorithm  upon  which  MELVIN  i  based  consists  of 
two  parts: 


T1  ires  ho  Id  De  t  e  nil  i  na  lion:  -  -  D  i  v  i  de  the  measured  PTVIR  into 
four  equal-duration  intervals.  For  each  of  the  four 
int  ervals  calculate  tin  average  power  using,  the  1  ir.st  J  T 
consecutive  samples  Select  the  lowest'  average  power  o  I 
the  four  and  set  a  threshold  level.  I  ivc  t  imes  (7  dB 
above)  this  noise-power  averag,e. 


Delay  Spread  Calculation: - -Determine  the  peak  signal 
power  of  the  measured  PTVIR.  Begin  the  calculation  of 
che  mom  .nts  {mk}  when  a  PTVIR  sample  reaches  one-quarter 
of  that  peak  for  first  time.  When  the  PTVIR  drops  below 
half  threshold  for  the  first  time,  begin  calculation  of 
che  10-sample  average  power,  but  re- initializing  the 
calculation  if  the  signal  rises  above  twice  threshold. 

Continue  calculating  the  moments  until  two  successive 
10-sample  averages  are  less  than  half  threshold. 

The  highly  complex  and  stochastic  character  of  the  experimental  PTVIRs 
has  led  to  several  algorithms  in  the  search  for  consistent  measures  of  delay 
spread.  For  example,  Cox  [13]  also  used  the  radius -of- gyration  measure  but 
rejected  all  responses  30  dB  or  more  below  the  peak.  The  SRI  [9]  delay- spread 
measure  is  based  on  the  time  interval  between  the  first  and  last  crossing  of  a 
threshold  four  (4)  times  above  the  average. 


The  delay-spread  calculations  are  summarized  in  Tables  3-3  (Vertical 
Polarization)  and  3-4  (Horizontal  Polarization)  and  grouped  according  to 
ante  na  height  (trunk  or  canopy),  season  (summer  or  autumn),  path  length  (927, 
471,  and  361  feet),  and  frequency  (451,  751,  1251  and  1751  MHz).  Column 
headings , abbreviated  for  compactness,  are  as  follows: 


Tp/Fi 

Date 

Range 

Ht 

Ant 


Tape  and  File  Numbers  (SRI  Archive) 

Date  of  data  acquistion 

Transmitter/Receiver  Antenna  Separation 

Transmitter/Receiver  Antenna  Height 

Transmitter/Receiver  Antenna  Type: 

L  —  log- periodic 

0  -  Omni-directional 

Upper  case  —  vertical  polarization 

Lower  case  —  horizontal  polarization 


Channel  Noise 
Nmin,  fsig 

Delay  Spread 
T 

ex;: 


Estimated  noi.se  power  floor  of  the  PTVIR 
(relative  to  peak  power)  and  its  .stan¬ 
dard  deviation  (relative  to  noise  floor) 

Delay  Spread  Moments  (Episode  Average) 
Excess  Delay  [ Eq .  (3-8)] 


-v" 


Rail  i  us  -  o  t  -  Cyra  t  i  on  [  Eq  .  (3-3) 

Do  1  ay  Sp read  [ Eq .  ( 3  -  9  )  ] 

Standard  Deviation  of  T 


Episode  averages  and  standard  deviations  are  based  on  seven  (7)  PTVIRs  taken 
from  successive  data  sets  (refer  to  Figure  3-  1.2)  and  separated  in  t  ime  by  16.  7 


ii!  i  1  1  i  seconds  . 


I 


i 


4 


20 


3 
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Table  3-2  : 


Delay- Spread  Calculations  (Vertical  Polarization) 


[Coventry,  Connecticut  -  1987] 


Tp/Fi 

Date 

Range 

(ft) 

Ht 

(ft) 

Ant 

Freq 

(MHz) 

Channel 

Nmin 

(dB) 

Noise 

?fl? 

T 

exc 

l - 

Delay  Spread 

T  T 

r  o  g  spr 

■nano'saconds- 

a 

--!? 

Trunks 

(Summer) 

T225F1 

09/01 
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Table  3-3: 


Delay- Spread  Calculations  (.Horizontal  Polarization) 
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Tp/Fi 

Date 

Range 

(ft) 

Ht 

(ft) 

Ant 

Freq 

(MHz) 

Channel 

Nmin 
(  dB  ) 

Noise 
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T 
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T 
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12 

11 

451 

-36.2 
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II 

11 

It 

II 
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24 

48 

72 

18 
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II 

II 

17 

II 
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40 

48 

88 

37 
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II 

II 

>1 

11 
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- 

- 

- 

- 
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20 

11 
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24 

24 

44 

0 
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II 

II 

II 

II 
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20 

32 

52 

4 
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II 

II 

II 
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44 
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52 
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36 

7*  U 

9 

T273F4 

11/10 

II 

It 

" 

699 

•38.8 

0 . 4 

12 

2  8 

40 

10 

T285F7 

n/1.1 

M 

” 

ll 

699 

-40.5 

-0.6 

8 

28 

36 
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tt 

1  751 

-25.9 

-0.2 

12 

52 

4  4 

13 

f '  86  F4 

11/11 

" 

” 

1“ 

1751 

-  36 . 2 

0 . 3 

20 

3  4 

64 

4 

3 


9 


4 . 0  Cone lus ions 


to 


9 


$ 


e 


# 


Before  attempting  to  draw  any  meaningful  conclusions  from  the  delay 
spread  calculations  of  Tables  3-2  and  3-3,  it  seems  prudent  to  consider  first 
the  quality  (or  significance)  of  the  data.  Two  measures  of  data  quality  are 
provided  in  the  Tables  by' :  the  estimated  channel  (background)  noise  power 
(N  )  and  the  standard  deviation  of  the  delay  spread  (a  )  .  A  cursory 

min  spr 

review  of  the  Tables  quickly  reveals  the  strong  correlation  between  these 
two  calculated  variables:  the:  higher  the  channel  noise,  the  larger  the 

standard  deviation.  It  is  also  apparent  that  the  least  reliable  data  are 
associated  with  vortical  polarization,  high  frequencies  and  long  ranges,  and 
so  correlate  directly  with  greater  transmission  losses  and  smaller  received 
signals  associated  with  these  propagation  conditions 

4 , 1  Trunk -Dominant  Measurements 

Tables  3-2  and  3-3  clearly  reveal  that  for  the  trunk- dominant  propagation 

the  delay  spread  is  greater  for  vertical  polarization  than  for  horizontal 

polarization,  although  the  difference  appears  to  decrease  with  increasing 

frequency.  The  frequency  dependence  of  the  delay  spread  depends  upon  the 

polar Lr. ation:  decreasing  with  increasing  frequency  for  vertical  polarization, 

and  increasing  wi'h  increasing  frequency  for  horizontal  polarization  For 

both  vertical  and  horizontal  polarization,  the  delay  spread  increases  with 

increasing  path  length.  This  is  also  apparent  from  Figure  4-1  which  provides 

a  comparison  with  South  Perry  data  acquired  previously  (8  Figure  4 -45] .  The 

differences  In  delay  spread  (greatest  between  400  and  451  MHz)  can  be 

attributed  to  the  different  biophysical  parameters  characterizing  the  sites 

[3.3  7  stcms/acre  at  South  Perry  versus  769  stems/acre  at  Coventry;  9.4  inches 

average  tree -trunk  diameter  at  South  Perry  versus  4.6  inches  at  Coventry] . 

The  near'!  y- equal  first  and  second  moments  of  the  delay  spread  [T  and  T  , 

<-■  *  <•  1.  o  8 

respectively  ]  revealed  by  Figures  4-2  and  4-3  coni  inns  the  ne.i  r- euponon  t:  i  al 
dependence  predicted  by  theory  [Equation  (211)] 


4 . 2  Canopy -Dominant  Measurements 
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Figure  4-2:  Delay  Spread  Plots 
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dependent  on  frequency,  if  at  all.  Distance  dependence  cannot  be  inferred 

from  the-,  limited  data  of  the  Tables.  Figure  4-4  clearly  reveals  the  delay- 

spread  function  to  be  non- exponential  since  the  second- central -moment  T  is 

rog 

several  times  larger  than  the  first -moment  T  .  Especially  surprising  is  the 
observation  that  for  canopy- dominant  propagation,  delay  spread  seems  to  be 
greater  in  winter  when  the  trees  are  tare  (no  foliage)  than  in  summer.  Rain, 
snow  and  sleet  seem  to  have  no  effect  on  delay  spread. 
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APPENDIX  A: 


Program  TP-6 


A 


(This  page  intentionally  left  blank] 


A 


f  tn77 
$MSEG  5 
$ FILES (1,1) 

program  tp6 (,9b) 
ema  hufi , bufq , power 
integer  trimlen 
charac cer*1600  string 
integer  bl(8192) ,b2(8192) 

integer  bufi(8192)  , bufq (11192)  ,bufd,  sorti(Sll)  ,  sortq(511) 
real  power(511 , 16) 
integer  ihard(28) 

C  EPISODE  SECTION 
C  INDEX  IS  CHANNEL 
INTEGERS 


LC0DE(2) , 
LBAND(2) , 
LCYCL( 2 ) , 
IPA(2) , 


CODE  LENGTH  255,511,1023,2047-0,1,2,3 
BAND  WIDTH  H , M , L , CW-0 ,1,2,3 
CYCLES/SAMPLE  1 , 2 , 4 , 8-0 , 1 , 2 , 3 
POWER  AMP  ASSIGNMENTS 


REALM 


IFREQ(2.)  , 
J  FREQ ( 2 ) 


!  RCVR  FREQ  MHz 
!  TX  FREQ  MHz 


INTEGERS 

IATYP ( 2 ) 


REALM 


IAL0C ( 2 ) 


ADIR(2) , 
AHGT ( 2 ) 


INDEX-1  FOR  RECEIVER,  -2  FOR  TRANSMITTER 

!  ANTENNA  C0DE:1~0MNI  2-VERTICAL 
3-HORIZONTAL  4-RH 
!  5-LH  (RCVR  ONLY) 

!  -1  ANTENNA  ON  MAST,  =2  ON  ROOF, 

-3  ON  TRIPOD 

!  ANTENNA  DIRECTION 
!  ANTENNA  HEIGHT 


EQUIVALENCE (IHARD(l) ,LCODE) 

EQUIVALENCE (IKARD( 3) , LB AND) 

EQUIVALENCE (I HARD ( 5) , I FREQ) 

EQUIVALENCE (IHARD( 9) , IPA) 

EQUIVALENCE ( I HARD ( l 1 ) , J  FREQ > 

EQUIVALENCE (I HARD ( 15) , IATYP) 

EQUIVALENCE (IHARD( 17) , ADIR) 

EQUIVALENCE ( IHARD( 21) , AHGT ) 

EQUT.VALENCE(  IHARD(  25  )  ,  LCYCL) 

EQUIVALENCE ( iHARD(27) , IALOC) 

THE  REMAINING  4  WORDS  ARE  SPARES 
equivalence  (string, b2) 
equivalence  (  IHARD , B2 ( 7169 ) ) 


call  igbuffbl ,8192) 


call  ehpar 


READ( 7 , iostat -  ios,  err— 900, 
Length  -  rr im Lent s f ring; 


‘nd-450)  L2 


OnLN ( 10 , FILL- ’ TV I R . DAT  * ) 


" , F7 . 2 , "  RECEIVER  FREQ?  - 
TRANSMITTER  FREQ 2  —  ”  ,  F  7 , 2  )  '  .) 


WRITE( 10.' ( "RECEIVER  FREQ1  - 
A  /  "TRANSMITTER  FREQ1  ,F7 .2, 


&  I FREQ (1) , I FREQ ( 2 ) ,JFREQ(1) ,JFREQ(2) 


WRITE ( 10, ' ("RX  ANTENNA  HEIGHT-" , F7 . 2 , "  R X  ANTENNA  DIRECTION-" , F7 . 2 , 
&/"T X  ANTENNA  HEIGHT-" , F7 . 2 , "  TX  ANTENNA  DIRECTION-" , F7 . 2) ' ) 

&  AHGT(l) , ADIR( 1 ) ,ADIR(2) , AHGT(2) 


IF 

(IHARD(l) 

.  eq. 

0) 

Write (10 , ‘ 

’("Code  length 

-  255")') 

IF 

(IHARD(l) 

.eq. 

1) 

Write ( L0 , ’ 

("Code  length 

-  511")') 

IF 

(IHARD(l) 

.eq. 

2) 

Write(10 , ' ( "Code  length 

-  1023") ' ) 

TF 

(IHARD(l) 

.  eq . 

3) 

Write ( 10 , 1 

'("Code  length 

-  2047") ' ) 

if 

( ihard( 3 ) 

.eq. 

0) 

write ( 10 , 1 

’("Bandwidth  - 

250  MHz") ' ) 

if 

(ihard(3) 

.eq. 

1) 

write(10 , 1 

1  ( "Bandwidth  — 

125  MHz") ' ) 

if 

( ihard( 3 ) 

.  eq. 

2) 

write ( 10 , 1 

’("Bandwidth  - 

50  MHz") ' ) 

if 

( ihard( 3 ) 

.  eq . 

3) 

write ( 10 ,'( "Bandwidth  - 

0  MHz") ' ) 

DO 

III  -1,2 

if  (iii  .eq.  1)  then 

if  (ihard(14+iii)  . eq .  1)  then 

write  (10, '("RX  antenna  type  -  Omni")') 
else 

if  ( ihard( 14+iii)  .eq.  2)  then 

write  (10, '("RX  antenna  type  -  Vertical")') 

else 

if  ( ihard(  1.4+iii)  .eq.  3)  then 

write  (10, '("RX  antenna  type  -  Horizontal")’) 
else 

if  ( ihard( 14+iii)  .eq.  4)  then 

write  (10. '("RX  antenna  type  -  RHC")') 

else 

if  ( ihard(14+iii)  .eq.  5)  then 
write  (10, '("RX  antenna  type  -  LHC")') 

end  if 
end  if 
end  if 
end  if 
end  if 
else 

if  ( ihard( 14+iii)  .eq.  1)  then 

write  (10, '("TX  antenna  type  -  Omni")') 
else 

if  ( ihard( 14+ii i)  . eq .  2)  then 

write  (10, '("TX  antenna  type  -  Vertical")') 

else 

if  ( ihard( 14+iii )  .eq.  3)  then 

write  (10, '("TX  antenna  type  -  Horizontal")') 
else 

if  ( ihard( 14+iii)  .eq.  4)  then 
write  (10, '("TX  antenna  type  -  PNC")') 

end  if 
end  if 
end  if 
e  nd  i  f 
END  IF 
END  DO 

DO  ircvr  ~  1,2  !  Do  for  both  receivers 

do  iset  -1,16  !  each  set  has  8  sweeps  that  will  be  averaged 

DO  10  J  — 1,4  !  j  -  channel,  no. 


A 


4 


(bufi(k) , k-1 , 8192 ) 
(bufq(k) ,k-l, 8192) 


if(  j  .eq.  1  )  Chen 

READ(7  .  iostat-ios  ,  err-900  ,  end— 950) 

READ ( 7 , ios  tat-ios , err-900 , end-950) 

iloc  -  6  *  1024 
DO  it  =  1,511 

sorti(it)  -  bufi(it  +  iloc) 
sortq(it)  -  bufq(it  +  iloc) 

End  do  !  end  it 
call  shell_sort(sorti , 511) 
call  shell_sort(sort.q ,  511) 
iavi  -  sorti(256) 
iavq  —  sortq(256) 

!  end  DC  Bias  determination 

do  m  —  1 ,  511 
x  - 

1  (( float (bufi (m+  iloc))-  float(iavi) ) )**2 .0  + 

2  (( float (bufq (m+  iloc))-  f loat( iavq) ) )**2 . C 

if(  x  . gt .  0.0  )  then 

power (m , iset)  -  10 . 0*ALOG10 (x) 
else 

power (m , iset)  -  -99. 
end  i  f 
end  do  !  m 

else 

read(7)  bufd 
read(7)  bufd 
end  if 
continue 

END  DO  !  iset  do  loop 
do  20  i-1,511. 

wr i te ( 10 ,  '  ( i5 , 16f 7 . 2 )  ' )  i , (power ( i , iset ) , I se t-1 , 16) 

END  DC  !  ircvr  do  loop 

1  STOP  -  LENGTH/ 80 
DO  IS  -  l.ISTOF 

write  ( 10  .  '  (a80/)  '  )  strings  (IS-  1 )  *80r-L  :  IS*80) 

END  DO 

I F( I St0P*80  . LT.  LENGTH)  wr i te ( 10 , ' ( a80/) ' ) 

■N  s tr ing ( i stop  •804 1 : length) 

c lose (10) 
s  t  o  p 

write) 1 , ' ( "Error  on  read  no .  — " , i4) '  )  ios 

stop 

write)  1,'("ECF  on  read")'  ) 
e  nti 


A 
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Program  TRflSFR 
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C2.34567 

G 

PROGRAM  TRNSFR 

C  THIS  PROGRAM  READS  CECOM  ASCII  FILES  FROM  TAPE,  PARTITIONS  THEM, 

AND  RE-WRITES  THEM  ON  DISK 


CHARACTER* 10 9  TXT 
DIMENSION  PWR (511, 16) 

5  FORMAT (109A.) 

10  FORMAT ( 5X, 16F7 . 2 ) 

11  FORMAT (8F7. 2) 

DO  100  1-1  TO  48 

READ (10, 5)  TXT 

WRITE (6, 5)  TXT  [Display  Header  on  Screen] 


100 

WHITE (11, 5)  TXT 

DO 

300  J-1,511 

300 

READ (10, 10)  (PWR(J.K), 

K-1,16) 

DO 

400  J-1,511 

400 

DO 

WRITE (12, 11)  ( PUR ( J , K ) 
401  J-1,511 

,  K-1,8) 

[  TVl’Rs  1-8,  RCVR  1 

401 

WRITE ( 14 , 11)  ( PWR ( J , K ) 

,  K-1,8) 

[ TVIRs  1-8,  RCVR  2 

STOP 

END 


/*  TRNSFR  EXEC  */ 

DO  I  -  l  TO  8 
J  -  2*1-1  ■ 

' FILEDEF  KEVIN1  TAPI  NL'  J 
' FILEDEF  SEVIN1  DISK  TEST  FILE  B' 


'MOVEFILE 

KEVIN1 

SEVIN1 ' 

[ Transfers 

file  from 

tape 

:  o  d  i  s  k 

' ASCTOEBC 

TEST  FI 

LE  B  TEST 

EBC 

B' 

[Converts  fil 

e  from  ASC 

II  to 

EBCDIC 

’ FILEDEF 

10  DISK 

TEST  EBC 

B' 

' FILEDEF 

li  DISK' 

•F'  |  |  I  !  | 

'  HDR ' 

'  DATA 

B  ( L.RECL  109  ' 

[ File  H 

eade  r 

' FILEDEF 

12  DISK’ 

'  F'  mu 

'A'  ' 

DATA 

B' 

[ PTVIRs 

1  -8  . 

RCVR  i. 

FILEDEF 

L  3  DISK' 

'  f'  mu 

'  B  '  ' 

DATA 

B' 

:  l’TVIRs 

“  •  1  6 

RCVR  1 

' FILEDEF 

13  DISK' 

'F'  1  UN 

DATA 

B’ 

[ PTVIRs 

1  -8  , 

RCVR  2 

*' FILEDEF 

14  DISK' 

'  F '  1  1  I  1  1 

' D '  ' DATA 

B' 

[ PTVIRs 

9  -  1  (> 

RCVR  2 

' FORTGO  TRNSFR' 

[  Ex  emit  es 

Fortran  Pr 

ogram 

TRNS FR 

'tape  r:-;w ' 


I 

i 
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APPENDIX  C: 

Program  ECK'JP 


10  !  BCKUP  ~  On  Blue  ?  -  3/15/90 

2  0  S 

30  OPTION  BASE  1 

40  DIM  Suf$ (3) [8] 

•50  DATA  "KDRA  » . ”AA  ",»CA 

60  READ  Suf $ ( *7 . 

70  FOR  I®1  TO  4 

SO  FOR  J=1  TO  3 

90  Sor$="F"&VAL!  'I)  &Suf  $  ( J) 

100  Sor$«Sor$  [3,10] 

•  110  r;<P.S$=JtT22  6F"  TVAL$  (I  MOD  9)&Suf$(J) 

120  Des$=Des$ [ 1 , 7 ] & "A" 

i-»0  PRINT  Sor$  ,  "  l!  ,  Des$ 

140  COP\  Sor$&” :HP8290X, 700,0"  TO  Des$&" : CS80 , 70 

150  NEXT  J 

160  PRINT 

•  170  NEXT  I 

180  END 
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D  -  - 


* 


a 


10 
20 
30 
4  0 
50 
60 
70 
80 
90 
100 
1.10 
12  0 
13  0 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
230 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
44  0 
450 

4  60 
470 
480 
490 
500 
510 
520 
530 
540 

5  50 
560 
570 
580 
5  90 
600 


!  KELVIN  -  On  Blue  #??  -  3/7/9C 

; 

!  This  Program  Reads  an  ITS  ASCII  Data  '  ile  From  Disc,  Rewrites  It 
!  as  a  STAT  Package-Compatible  BDAT  Data  File  on  disc,  and  produces 
!  Summary  Data  Sheets.  [Delay  Spread  is  not  calculated.] 

i 

OPTION  BASE  1 
REAL  Sc (20) 

DIM  T$ [80] , Vn$ (50) [10] ,Sn$(20) [10] 

DIM  A$(50) [74] ,B$(513) [74] ,Y(8,511) ,Tp$(6) [4] ,Rc$(4) [2] 

DATA  "T247" , "T226" , "T225" , "T273" , "T264" , "T285" 

DATA  "Rl" , "R3 " , "R2" , "R4" 

READ  Tp$  (■■’•■)  ,  Rc$  (  * ) 

C=1 

INPUT  "ENTER  the  number  of  the  file  to  be  copied", C 
Cc=C  MOD  3 
Tc=C  DIV  8+1 
PRINT  CHR$ ( 12 ) 

PRINT 

PRINT 

MASS  STORAGE  IS  " : HP82901 , 700 , 0" 

!  Reading  the  Header 

C$— " F" &VAL$ (C) & "HDRA _ " 

ASSIGN  @Inpath  TO  C$ 

ENTER  @Inpath;A$(*) 

GOSUB  610 

ASSIGN  @Inpath  TO  * 

INPUT  "Do  you  want  to  store  this  data  on  disc?  [Yes=l  No=0]",Yn 
IF  Yn=0  THEN  1410 
PRINT  CHR$ ( 12 ) 

i 

FOR  L=1  TO  4  STEP  2  !  Set  STEP  to  1  to  get  all  16  TVIRs 

C$=" F" &VAL$ (C) &CHR$ (64+L) &"A  " 

C$=C$ [ 1, 10] 

F$=Tp$ (Tc) & "  F" &VAL$ (Cc) &Rc$ ( L) 

DISP  "Reading  ASCII  Data  File  "&C$ 

ASSIGN  @ Inpath  TO  C$ 

ENTER  @Inpath;B$ (*) 

DISP  "Processing  ASCII  Data  from  File  "&C$ 

FOR  J=3  TO  513 
FOR  K=1  TO  Nv 

Y  (  K .  J-2) =10. 0A  ( VAL  (  B$  (J)  [7*K-4 , 7*K+2  ]  )/-l  0) 

NEXT  K 
NEXT  J 

ASSIGN  @Inpath  TO  * 

!  Write  First  Record  to  Disc 
DISP  "Writing  Header  on  Disc" 

MASS  STORAGE  IS  " : HP8290X, 700 , 1" 

CREATE  BDAT  F$ , INT ( ( 8 *Nv*No) / 128 0 ) +3 , 12 8 0 
ASSIGN  @F ilel  TO  F$ 

OUTPUT  @ Pilel  ;T$  ,  No ,  ilv,  Vn$  ( * )  ,  Ns,Sn$(*)  ,Sc(*) 

!  Write  TVIRs  (1-8)  on  Disc 

DISP  "Writing  PTVIRs  to  Disc  -  "&F$ 

OUTPUT  @Filei , 2 
OUTPUT  @Filel;Y(*) 

ASSIGN  @Filel  TO  * 

MASS  STORAGE  IS  " : HP8290X, 700 , 0" 

DISP  "  " 

NEXT  L 
GOTO  1410 


D 


« 


610  !  Constructing  the  Header 

620  DATA  " PTVIR1 " , " PTVIR2 " , " PTVIR3 " , " PTVIR4 " , "PTVIR5" , " PTVIR6 " 

630  DATA  "PTVIR7", "PTVIR8", "  " , " ? " , " ? " 

640  FOR  1=1  TO  15 
650  READ  Vn$ (I) 

660  NEXT  I 

670  Vn$  (12)  =A$  (3)  [60, 65] &"  87"  S  Year  must  be  changed  'manually 

680  Vn$ (16) =VAL$ (L) 

690  Vn$ (17) =TRIM$ (A$ (47) [POS (A$ (47) , "Code  length  =  ")+14;5]) 

700  Vn$  ( 18 )  “”TRIM$  ( A$  (48)  [  POS  (A$  (48 )  ,  "Bandwidth  =  ")+14;4]) 

710  Vn$ (19) =TRIM$ (A$ (43) [POS(A$(43) , "FREQ1  =")+7;7] ) 

720  Vn$ (20) =TRIM$ (A$ (50) [POS (A$ (50) , "TX  antenna  type  =  ")+19;l]) 

730  Vn$ (21) ="  " 

740  Vn$ (22) =TRIM$ ( A$ ( 4 6 )  [ POS ( A$ ( 46 )  , "TX  ANTENNA  DIRECTION*" ) f 2 1 ; 7 ] ) 

750  Vn$ (23 ) =TRIM$ (A$ (46) [POS (AS (46) , "TX  ANTENNA  HEIGHT=" ) +18 ; 7 ] ) 

7  6C  Vn$  (24)  =TRIM$  (A$  (49)  [POS  (A$  (49)  ,  "RX  antenna  type  =  ")+19;l]) 

770  Vn$ (25) " 

780  Vn$  (26)  =TRIM$  ( A$  ( 4  5)  [  POS  ( A$  ( 4  5 )  ,  "RX  ANTENNA  DIRECTION="  )  4-2  1 ;  7  ]  ) 

790  Vn$ (27) =TRIM$ (A$ (45) [POS (A$ (45) , "RX  ANTENNA  HEIGHT*" ) +16 j 7 ] ) 

300  Bw»VAL(Vn$ (18) ) 

810  Nv=8 

820  No=VAL(Vn$ (17) ) 

830  Ns=0 

840  Fl$=Tp$ (Tc) &"F"&VAL$ (Cc) 

850  T$="Tape&File:  "&F1$&"  Date:  "&Vn$(12)&"  Time:  "&Vn$(13)&"  Experiment. 
"&Vn$ ( 14 ) &" ( "&Vn$ ( 15) &") " 

860  DISP  "" 

370  MASS  STORAGE  IS  " : HP8290X, 700 , 0" 

880  Printo=0 

890  PRINT  TAB (31), "WPMS  File  Description" 

900  PRINT 

910  PRINT  TAB (7) , T$ 

920  PRINT 
930  PRINT 
940  PRINT 

950  PRINT  TAB (10) , "Channel:  "£Vn$(16) 

960  PRINT 

970  PRINT  TAB(10) , "Code  Length:  " &Vn$ ( 17 ) &"  Chips" 

980  PRINT 

990  PRINT  TAB (10) , "Bandwidth:  "&Vn$(18)&"  MHz" 

1000  PRINT 

1010  PRINT  TAB(10) , "Carrier  Frequency:  "&Vn$(19)&"  MHz" 

1020  PRINT 

1030  PRINT  TAB ( 10) , "Transmitting  Antenna:" 

1040  PRINT 

1050  PRINT  TAB ( 18 ), "Type:  "&Vn$(20) 

1060  PRINT  TAB (18) , "Location:  "&Vn$(21) 

1070  PRINT  TAB (18) , "Direction:  "&Vn$ (22 ) &"  Degrees" 

1080  PRINT  TAB (18) , "Height:  "&Vn$(23)&"  ft" 

1090  PRINT 

1100  PRINT  TAB( 10) , "Receiving  Antenna:" 

1110  PRINT 

1120  PRINT  TAB (18) , "Type:  "&Vn$(24) 

1130  PRINT  TAB( 18) , "Location:  "&Vn$(25) 

1140  PRINT  TAB ( 18 ), "Direct ion :  "&Vn$(26)&"  Degrees" 

1150  PRINT  TAB ( 18) , "Height:  "&Vn$(27)&"  ft" 

1160  PRINT 

1170  PRINT  TAB ( 10) , "Path  Loss:" 

1180  PRINT 

1190  PRINT  TAB ( 18) , "Derived  from  WRSL:  "&Vn$(30)&"  dBm" 


€ 


4 


u 


4 


1200  PRINT  TAB ( 18 ), "Derived  from  PTVIR;  "&Vn$(31)&"  dBm" 

1210  PRINT 

1220  PRINT  TAB(10) , "Received  Power:" 

1230  PRINT 

1240  PRINT  TAB (18) , "Wideband:  "&Vn$(32)&"  dBm" 

•  1250  PRINT  TAB ( 18 ), "Narrowband:  "&Vn$(33)&"  dBm" 

1260  PRINT  TAB (18) , "PTVIR:  "&Vn$(34)&"  dBm" 

1270  PRINT 

1280  PRINT  TAB ( 10 ), "Transmitted  Power  (EIRP) :  "&Vn$(35)&»  dBm" 

1290  PRINT 

1300  IF  (Printo=l)  THEN  GOTO  Final 

£  1310  INPUT  "Do  You  Want  This  Information  on  Hard  Copy?  [Yes=l  No=0]",Pyn 
1320  IF  (Pyn~0)  THEN  GOTO  Final 
1330  PRINTER  IS  9 
1340  PRINT  CHR$ ( 13 ) 

1350  Printo=l 
1360  GOTO  890 

0  1370  Final: IF  (Printo=0)  THEN  GOTO  1400 
1380  PRINT  CHR$ ( 12 ) 

1390  PRINTER  IS  1 
1400  RETURN 
1410  END 


D 
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MELVIN  on  Blue  #??  -  3/12/90 


m 


10 
2  0 
30 
40 
50 
6  0 
70 
80 
90 
100 
110 
120 
13  0 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
24  0 
250 
260 
27  0 
280 
290 
300 


m 


310 

320 

33  0 

34  0 

3  50 
360 
370 
33  0 
.390 

4  00 
410 
4  2  0 
43  0 
440 
4  50 
4  60 
■17  0 
430 

4  9  0 
500 
510 
520 

5  3  0 
5  4  0 
c.  5  0 
5  6  0 
57  0 
5  S  0 


!  This  Program  Reads  SRI  BOAT  Data  Files  from  Disc,  Calculates 
i  the.  Delay  Spreads  (second  central,  moments)  using  the  Altman 
!  Algorithm  (CyberCom  Technical  Report  CTR-117-01,  Section  4.3.3), 

!  and  produces  Summary  Data  Sheets  (op.  cit..  Table  4-4). 

j 

OPTION  BASE  1 

INTEGER  PI , Pp ( 4 ) , Dp , P , Q , Rf 
REAL  Sc (20) 

DIM  T$ [ 80  J , Vn$ ( 50 )  [10] , Sn$(20)  [10] ,Nn(20)  ,Titlel$[30]  ,Title2$ [ 80 ] 
F$="T225F6R1" 

ASSIGN  @Path  TO  F$ 

ENTER  @Path ;T$ , No , Nv, Vn$ ( *) ,Ns,Sn$(*) ,Sc(*) 

Cf=1000/VAL(Vn$ (13) )  !  Sample  interval  in  nanoseconds 

Nc=C 

i  Noise  Window  Definition 
DATA  12,40,5 
READ  PI , Dp , Fac 

i  INPUT  "Noise-Window  Start  and  Width  (in  Samples) ?" , PI , Dp 
!  INPUT  "Threshold/Noise  =  ?  [Default  -  5]", Fac 
Pp(l)=Pl  !  Location  of  fi.rst  noise  window 

FOR  J=2  TO  4 

Pp(J)=Pp(J-i) +lNT(No/4)  !  Location  of  subsequent  noise  windows 

NEXT  J 

i 

ENTER  @ Path, 2 
ALLOCATE  Z (No) 

ALLOCATE  Pkl (Nv) , Avq (Nv) , Sig (Nv) , Sprd (Nv) 

PRINT  TAB (11) , "PTVIR  Nmin  Nsig  Start  Mean  Sigma  Spread 

FOR  Typ-1  TO  Nv  !  PTVIR  loop 

!  Read  PTVIRs  from  Disc 
ENTER  @Path;Z (*) 

MAT  SEARCH  Z(*),LOC  MAX;Pk 
Zpk~Z (Pk) 

Inc-INT  (  (No+l'i/S) 

Rf-Pk-Pk  MOD  Inc+l-Inc 
IF  F.f=l  THEN  44  0 

ALLOCATE  Zc(No)  !  Shifts  PTVIR 

MAT  Zt  =  Z 

MAT  Z ( 1 : No—Rf +1 ) =  Zt(Rf:No) 

MAT  Z (No  -Rf +  2 : Mo) =  Z t ( 1 : R r - 1 ) 

DEALLOCATE  Zt(*) 

!  Finding  rms  noise  within  the  windows 
DATA  0,0,0,1,1,100000 
RESTORE  *150 

READ  SO , SI , S2 , Pk , Flag , Nmin 
FOR  P-1  TO  No 

IF  Z(Pk)<Z(P)  THEN  Pk=P  !  Locates  location  of  PTVIR  max. 

IT  FLat,=  5  THEN  67  0 

IF  P<Pp (Flag)  THEN  670  !  Rejects  non-window  noise  samples 

SO-SCfl 

si«si+Z(F) 

S2-S2+Z (P) *Z ( P) 

IF  P>Pp (Flag) +Dp  THEN 
N ( Flag; =INT(S1/S0) 

Mm=S0*S2-Sl*Sl 
T  c  N  ( F 1  a  g )  <  N  m  i.  n  T  HEN 
Ni=Flag 

B  -  3 


< 


600 

Nmin=N (Flag) 

610 

Nsig-INT (SQR (Mm) / 

620 

END  IF 

630 

RESTORE  450 

640 

READ  SO , SI f  S2 

650 

Flag=Flag+l 

660 

END  IF 

670 

NEXT  P 

680 

GOSUB  2030 

690 

GOTO  720 

700 

INPUT  "Retain  for  Delay 

710 

IF  Yn~0  TEEN  1330 

72  C 

Nc=Nc+l 

730 

!  Calculation  of  Delay 

740 

RESTORE  450 

7  50 

READ  S 1 , S  2 , S  3 

7  60 

RESTORE  450 

770 

READ  Fig, Fig, Fcg 

78  0 

S4  =  0 

7  90 

Jj=l 

800 

P2  =  l 

810 

Ct.  =  10 

820 

Th=Fcic*Nmin 

830 

FOR  T-.l  TO  10 

840 

Nn ( I) -0 

850 

NEXT  1 

8  6  0 

FOR  Q=1  TO  No 

370 

P=Q 

880 

IF  Fog-1  THEN  1.2 5 G 

890 

IF  Flg=l  THEN  960 

900 

11  Z(P)>Zpk/4  THEN 

910 

Flg=l 

920 

Pk3 (Typ)=P 

930 

ELSE 

940 

GOTO  1260 

950 

END  IF 

960 

IF  Fig=l  THEN  990 

Q70 

IF  Z(P) < . 5*Th  THEN 

980 

Fig-l 

990 

IF  Z ( P) >2 *Th  THEN 

1000 

Jj  =  l 

1010 

Fig=0 

1020 

Su=0 

10  3  0 

P2=-U 

104  0 

GOTO  1200 

1050 

END  IF 

10  6  0 

P2-P2+1 

107  0 

Su=Su+Z { P) 

1080 

IF  P2  >Ct  THEN 

109  0 

Nr (J  j ) =Su/Ct 

1100 

IF  J  j  -•=  1  THEN  1160 

i  :i  .1  o 

IF  Nn (J j ) < . 5*Th  A 

112  0 

Pend=P 

1 1  3  0 

Fog  - 1 

114  0 

GOTO  1240 

1150 

END  IF 

116  0 

117  0 

S  a-  0 

1 1. 3  0 

P  2  ~  0 

1190 

END  IF 

Re-initialize  moments 
Identifv  next  window 


Calculation?  [Yes=l  No=0] 


Yn 


Spread 


■a  arts  Integration 


Nn  { 3  j  ~  .1 )  •  .  5  *Th  THEN 


( 


1 


i 


I 


4 


4 


« 


<1 


0 


•i 


9 

1200  END  II 

1210  S1=S1-+I(P) 

1220  S2=S2+Z  ( P)  *Q 

1230  53=S3+Z (P) *Q*Q 

1240  IF  Z (P) <Fac*Th  AND  Fog=l  THEN  1260 

#1250  S4=S4+Z ( P) 

1260  NEXT  Q 

1270  Avg  (Typ)  — S2/S1 

1280  IF  Pk<Mid  THEN  1300 

1290  Avg  (Typ) -Avg  (Typ) -f-M id 

13  00  Sig (Typ) =3QR (S3  *S1-S2  *S2 ) /SI 

#1310  Sprd(Typ) =INT (Avg (Typ) -Pkl (Tvp) +Si g(Typ) ) 

1320  PRINT  USING  Frmtl ;Typ , Nrcin , Nsig, Pkl (Typ) , Avg (Typ) , S ig (Typ ) , Sprd (Typ) 
1330  NEXT  Typ  i  Next  PTVIR 

1340  Frmtl : IMAGE  11X, DD, 2 (4X, DDDDD) , 4 ( 5X, DDDD) 

1350  Strt=f UM ( Pkl ) /Nc 
1360  Dspr=SUM(Sprd)/Nc 
#1370  PRINT  USING  Frmt2 ; Strt , Dspr 

1380  Frmt2: IMAGE  8X, "Averages" , 2 OX, DDDD, 2 3X, DDDD 
1390  WAIT  2.5 
1400  ! 

1410  ASSIGN  @Path  TO  * 

1420  DISP  "" 

^  1430  GRAPHICS  OFF 
1440  GOSUB  1540 

1450  INPUT  "Do  you  want  a  hard  copy?  [Yes-1  or  No=0]",Yn 
1460  PRINT  CHR$ (12) 

1470  GRAPHICS  ON 
1480  IF  Yn- 0  THEN  2220 
1490  PRINTER  IS  9 
1500  GOSUB  1540 
3.510  PRINT  CHR$  (12) 

1520  PRINTER  IS  1 


153  0 

GOTO 

2220 

154  0 

!  Pri; 

nt  File 

Description 

m  ]  550 

PRINT 

TAB (3 11 

, "WPMS  File  Description" 

9  1560 

PRINT 

1570 

PRINT 

TAB (10) 

,  T  i  1. 1  e  2  $ 

1580 

PRINT 

1590 

PRINT 

TAB (10) 

,  "  Chan  tie  1.  :  "  &  F$  [8,8  j 

1 6  G  0 

PRINT 

-  1610 

IQ 

PRI  vil 

TAB  ( 10 ) 

, "Code  Length:  "&Vn$(17)&"  Chips" 

162  0 

PR  INI 

L  6  3  0 

PRINT 

T  '/•  B  ,  1  o) 

,  "  Bandwidth  :  "iVn$  (  18) &"  Mil.:" 

164  0 

PRINT 

1650 

PRINT- 

TAB (10) 

, "Carrier  Frequency:  "&Vn$(19)&"  MHz" 

.1660 

PRINT 

m  167  0 

PRINT 

TAB (10) 

, "Transmitting  Antenna: " 

9  1680 

PRINT 

169  0 

PRINT 

TAB (18) 

,  "Type:  "SP/nS  (20) 

17  00 

PRINT 

TAP ( 18) 

.  !*  Locat.  ion  :  '■  S  Vn$  (  2  !  ) 

1.710 

PRINT 

TAB ( 18) 

, "Direct  ion :  "&Vn$ (  22  ) 

i&"  Degrees" 

1  7  2  0 
1.7  3  0 

PI!  I  NT 
PRINT 

TAB ( 18) 

, "Height:  "&  i$(23)&" 

f  t " 

1  7**0 
1750 

PRINT 

PRINT 

'1  AB  (  1  0  ) 

, "Receiving  Antenna: " 

1760 

PRINT 

TAB  (  1 8  ) 

, "T^pe :  " & V  n  $ ( 2  4 ) 

1  7  ?  0 

PRINT 

TAB ( 18) 

,  "  Locat  l  or  :  "  &Vn$  (  2  5  ) 

1  7  8  0 

PR  TNT 

TAB ( 18 ) 

,  ' 1 D  i  r  e  <:  t  v  o  n :  "&VnS  (  7  6 1 

i&"  !Vn|i  ees" 

1790 

p.R  if  r 

TAB (18) 

, "height:  •’ &Vn$  (  2  7  )  &  " 

1 1" 

gmufl 


1800  PRINT 

1810  PRINT  TAB ( 10 ), "Path  Loss:" 

1820  PRINT 

1830  PRINT  TAB ( 18 ), "Derived  from  WRSL:  "&Vn$(30)&"  dBm" 

1840  PRINT  TAB ( 18 ), "Derived  from  PTVIP:  "&Vn$(31)&"  dBm" 

1850  PRINT 

1860  PRINT  TAB ( 10 ), "Received  Power:" 

1870  PRINT 

1880  PRINT  TAB ( 18 }, "Wideband:  "&Vn$(32)&"  dSm" 

1890  PRINT  TAR (18) , "Narrowband:  "&Vn$(33)&"  dBm" 

1900  PRINT  TAB (18) , "PTVIR:  ”&Vn$(34)&"  dBm" 

1.910  PRINT 

1920  PRINT  TAB ( 10) , "Transmitted  Power  (EIRP) :  " &Vn$ ( 3 5 ) & "  dBm" 

1930  PRINT 

1940  PRINT  TAB( 10) , "Delay-Spread  Calculations:" 

1950  PRINT 
1960  Cf=l 

1970  ! PRINT  USING  Frmt4 ? Avg*Cf , Dspr*Cf 
1980  PRINT  USING  Frmt5 ;3trt*Cf , Dspr*Cf 

1990  FrmtA; IMAGE  17X, "Delay:  Average  =  "f4D,"  Std.Dev.  =  ",4D,"  ns" 

2000  Frmt5: IMAGE  24X, "Minimum  =  ",40,"  Spread  =  " ,4D,"  ns" 

2010  PRINT 
2020  RETURN 

2030  i  Plotting  of  Delay  Spread  Function 
2040  Frst=l 
2050  Lst-No+1 

2060  Titlel$="  "  !  "Delay-Spread  Function" 

2070  Title2$=T$[l,POS(T$, "Date") +14 ] 

2080  Crt-0 

2090  CALL  Plotpnts ( Z ( * ) , Frst , Lst ,Titlel$,Title2$,Cf,Crt) 

2100  GOTO  2210 
2110  ! GRAPHICS  OFF 

2120  INPUT  "Do  you  want  to  use  the  plotter?  [Yes=l,  No=0]",Ans 
2130  IF  Ans=0  THEN  2170 
2140  Crt=l. 

2150  INPUT  "When  the  plotter  is  ready  press  ENTER", Ans 
2160  GOTO  2  0Q0 

2170  !  INPUT  "Do  you  want  to  plot  a  sub-set.  or  quit?  [Sut-set-1,  Quit=0  ]  "  ,  Ans 

2180  ! IF  Ans=  0  THEN  1530 

2190  !  INPUT  "Enter  the  first  and  last  points  to  be  plotted.;  ",  Frst,  Lst 

2200  IGOTO  2930 
2210  RETURN 
2220  END 
2  2  3  0  ! 

2240  SUB  Plot  prts (array (*), First  pnt , Last  pnt , Main  _title$, Ref  t itle$ ,  Samp  int , 
Crt , OPTIONAL  01) 

2  2  5  0  ! 

2260  !  Array (*)  CONTAINS  THE  POINTS  TO  BE  PLOTTED 

2  2  7  0  ! 

2280  DIM  Dummy $ [ 80 ] 

2290  !  DISP  "Graph  Title:  [80  Char  Max,  DEFAULT**" &Ma  in  title$&" 

2  3  00  !  i_.IN  PUT  Dummy  $ 

2310  ! IF  Gummy $<>" "  THEN  Main  tit le$=Dummy$ 

2320  Labql : Answer$«"B" 

2330  !  INPUT  "Bar  (B)  or  Line  (  L)  Graph  ( DEFAULT-- B)  ?",  Answer? 

2340  'IF  Answer$<>" B"  AND  Answer$o"L"  THEN  Labql 
2  3 5 o  Smin* F i rst _pnt 

2360  Sola x= La st  pnt 

2370  ! DISP  "Plot  samples  between  Smin,Smax  [ DEFAULTS- ("; First  pnt ; " , " ; Last  pnt 

;  "  )  J " ; G CLEAR 


% 


# 


* 


f 


E  ~  6 


t'J  fJ  WfJIOWfJNJtJNJWKJMWWWrJWWtOMSJtOWMMWMWWfJIO  (V  t-J  M 


2380 
2390 
2400 
2  410 
2420 

•  2430 
2440 
2450 
2460 
2470 
2480 

•  2490 
2500 
2510 
2520 
2  53  0 
2540 

•  2550 
2560 
2570 
2580 
2590 
2600 

610 
620 
63  0 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
84  0 
850 
860 
370 
880 
890 
900 
910 
920 
930 
94  0 
9  50 


2  9  6  0 
•  2970 


! INPUT  Smin , Smax 
! IF  Smax<Smin  THEN  3110 
! First_pnt~Smin 
!  Last.  _pnt~Smax 
GINIT 

IF  Crt  THEN  PLOTTER  IS  703,"HPGL" 

GRAPHICS  ON 

Xr--.7  !  Reduced  from  .8 

Yr=.5  !  Reduced  from  .8 

cr=Xr 

IF  Xr>l  THEN  Cr=l 
CSIZE  4*Cr , . 6 
LORG  5 
LDIR  0 
PEN  1 

ON  ERROR  GOTO  Skipl 

VIEWPORT  18*Xr+15,  13 l*Xr+15 , 9 l*Yr+10 , 92 * Yr+10 !  SET  UP  AXES 
WINDOW  0,8,0, 1 
AXES  1,0, 0,0 

VIEWPORT  18*Xr+15, 131*Xr+15, 2 5* Yr+10, 26* Yr+10 
WINDOW  0 , 8 , 1 , 0 
AXES  1,0, 0,0 

VIEWPORT  16*Xr+15, 17*Xr+15,27*Yr+10,90*Yr+10 
WINDOW  1, 0, 0,  1 
AXES  0 ,  . 1 , 0 , 0 

VIEWPORT  132  *Xr+15 , 133*Xr+15, 27*Y.r+10 , 90*Yr+10 
WINDOW  0 , 1 , 0 , 1 
AXES  0, . 1, 0, 0 

VIEWPORT  4*Xr+15, 15*Xr+15 , 2 6*Yr+10 , 9 1* Yr+10 
WINDOW  0,1,-. 02,1.02 
CLIP  OFF 

FOR  1=1  TO  0  STEP  -.2 
MOVE  .5,1 

LABEL  USING  Imagel;! 

Imagel: IMAGE  MD.D 
NEXT  I 

VIEWPORT  18*Xr+15, 131*Xr-rl5, 9 ~*Yr+10, 97*Yr+10 
WINDOW  0 , 1 , 0 , 1 
MOVE  0 ,  .  5 
LORG  2 

LABEL  Ref  tit.le$  [  1 , 19  ] 

LORG  0 
MOVE  1,-5 

LABEL  Ref_title$ [20,34] 

VIEWPORT  0*Xr+15 , 4 *Xr+ 15 , 2 6 * Yr+10 , 9 1 * Yr+ 1 0 

WINDOW  0, 1 , 0,  1 

CLIP  OFF 

LDIR  PI/ 2 

MOVE  . 5 , . 5 

LORG  5 

LABEL  "Normalized  Power" 

LDIR  0 

VIEWPORT  18*Xr+15, 131*Xr+ 15 . 20*Yr+l0, 24*Yr+10 
WINDOW  0,1,0, 1 
MOVE  0 , . 5 
LORG  2 

LABEL  VAL$( First  pnt) 

MOVE  1 . , . 5 
LORG  8 

LABEL  VAL$ ( Last  pnt) 


2930  MOVE  .5,. 5 

2990  LORG  5 

3  000  .LABEL  "Sampler “ 

3010  Image2 : IMAGE  DDPDO 
3  02  0  LORG  2 

3030  VIEWPORT  18*Xrtl5, 13 l«Xr+15 , 10*Yr+10, 16*Yr+10 

3040  WINDOW  0, 1,0,1 

3050  MOVE  0 , . 5 

3060  UnitO=,,nsec" 

3070  LABEL  USING  Image4 ? Samp_int , Unit? 

3  08  0  Image4: IMAGE  "  Sample  Interval:  "  ,  MD.  3DE ,  IX,  R 

3  09  0  Yma:e»0  !  DETERMINE  Ymax 

3100  PSNUP 

3 11C  FOR  I=First  put  TO  Last  pnt-1 

3120  IF  ABS (Array (I) ) >Ymax “THEN  Ymax=AES (Array (I) ) 

3130  NEXT  I 

3140  VIEWPORT  18*Xr+15 , 131*Xr+15 , 14*Yr+10 , 20*Yr+10 

3150  WINDOW  0,1, 0,1 

3160  MOVE  0 ,  „  5 

3170  LABEL  USING  Image!; Ymax 

3180  Image 3 : IMAGE  "  Maximum  Power:  " , MD. 3DESZZ , "  (=1)" 

3190  VIEWPORT  18 *Xr+15 , 13  l*Xr+ 15 , 2 7  * Yr+10 , 90* Yr-f  10 

3200  IF  YmaxoO  THEN  WINDOW  First_pnt ,  Last_pnt ,  0 ,  Ymax 

3210  IF  Ymax=0  THEN  WINDOW  First_pnt , Last_ pnt , -1 , 1 

3220  !  PLOT  POINTS 

3230  MOVE  First_pnt,0 

3240  DRAW  Last_pnr,0 

3  250  IF  Answer$="B''  THEN  MOVE  Smin ,  0 

3260  IF  Answer$="L"  THEN  MOVE  Smin , Array ( Smin) 

3  27  0  FOR  I=Sm.in  TO  Smax-1 

3280  IF  Answer$="B"  THEN 

3290  MOVE  1,0 

3300  IF  A3S (Array ( I ) /Ymax) <.01  THEN  GOTO  3350 

3310  DRAW  I, Array (I) 

3320  ELSE 

3330  DRAW  I , Array ( I ) 

3340  END  IF 

3350  NEXT  I 

3360  IF  Smin-First  prt  THEN  3400 

3370  LINE  TYPE  5,2“ 

3380  MOVE  Smin, -Ymax 

3390  DRAW  Smin, + Ymax 

3400  IF  Smax=Last  pnt  THEN  3440 

3410  LINF  TYPE  5,2 

3420  MOVE  Smax, -Ymax 

34  3  0  DRAW  Smax ,  lYmux 

3440  LINE  TYPE  1 

3450  VIEWPORT  3*Xr+15, 146*Xr-M  5 ,  1 04  * Yr+10 , 110*Yr+10 

3  4  60  WINDOW  0,1, 0,1 

3470  CS1ZE  5*Cr,.6 

3 480  LORG  5 

3  4  90  MOVE  .5,.  5 

3500  LABEL  Main  title? 

3 510  X-LEN ( TR IM$ ( Ma i n  title?)) / 2 

3520  MOVE  . 5- (  6«5*Cr) / ( 14 1+Xr) *X, 0 

3 5 3 0  ! DRAW  .5+(.6*5*C r ) / ( 1  A  1  * X r )  *X, 0 

3540  Ski pi : OFF  ERROR 
3550  PENUP 

3560  PEN  0 

3570  ON  ERROR  GOTO  Subend 


m 


# 


# 


# 


8 


#5  SO  OX-Yxnax 
3590  Subend: SUBEND 


